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ABSTRACT 

Recent observations of tidal debris around galaxies have revealed that the structural 
properties of the spheroidal components of tidally disturbed galaxies are similar to 
those found in non-interacting early-type galaxies(ETGs), likely due to minor merging 
events that do not strongly affect the bulge region or to major mergers that happened 
a long time ago. We show that independently of merger events, tidal features like 
shells or rings can also arise if the the dark matter is an ultra light scalar field of 
mass ^ 10 _22 eV/c 1 2 * * * * * . In the scalar held dark matter (SFDM) model the small mass 
precludes halo formation below ^ 1O 8 M 0 reducing the number of small galaxies today, 
it produces shallow density profiles due to the uncertainty principle in contrast to the 
steep profiles found in the standard cold dark matter (CDM) paradigm, in addition to 
the usual soliton solution there exists dark matter haloes in multistates, characterized 
by ripples in their density profiles, which are stable provided that most of the halo mass 
resides in the ground state. We use the hydrodynamics code ZEUS to track the gas 
evolution in a background potential given by a superposition of the ground and first 
excited state of the scalar held, we study this configuration when it is initially unstable 
(excited state more massive than ground state) but by a population inversion in the 
states it eventually becomes stable, this could happen when haloes decoupled from the 
expansion of the universe and collapse to reach a state of equilibrium. We found that 
tidal structures like rings are formed at a particular radius as a direct consequence of 
the wavelike structure of the dark matter halo, thus they possess a certain degree of 
symmetry that can be used to distinguish them from remnants formed through the 
usual galaxy merging mechanism. In our simulations most of the gas concentrates in 
the center of the halo forming a dense and compact structure likely enhancing star 
formation in this region, however, some of the gas always remains trapped between 
the overdensity regions of the background SFDM halo predicting surrounding features 
even in non-interacting galaxies. In galaxies that are not strongly interacting with 
their neighbors, these structures may leave traces such as radial gradients in stellar 
metallicity, rings of dust, or even in the form of a younger population of stars at 
particular radii in contrast to the old stellar population dominating the bulge region, 
our results offer a new mechanism in which features like rings and shells can emerge 
purely from the quantum properties of this ultra light scalar held. 

Key words: circumstellar matter - infrared: stars. 


1 INTRODUCTION 

The standard model of cosmology (or CDM) assumes that 

the dark matter(DM) is cold and effectively collisionless. It 

assumes a hierarchical framework for galaxy formation and 

throughout their evolution galaxies are subject to frequent 

collisions and interactions with nearby galaxies that leave 

surrounding material in the form of tidal debris, therefore 
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studying tidal features can yield important clues on galaxy 
evolution. 


Tidal features can appear in the form of complicated 
structures ([Barnes fe Hernauistl Il992l h but it is frequent 
to observe struc tures like tidal t ails around the b ulge of 


early-type galaxies present shellsdMalin & Carter] 

1980, 

1983: iTurnbull.Brid 

e:es & Carterl 19991: van Dokkum 

2005; 

Schweizer & Seitzer 

19921: Michard & Prugniel 2004h which 


can be seen as smooth arcs of ellipses centered around 
its host galaxy decreasing their surface brightness for 
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larger r adii, there are also debris in the form of stellar 
streams (|Martmez-Delgado et al.ll2008 , 20091 , l201(ih . 

Several analyses of the spheroidal com pone nts of ti dallv 
disturbed early-type galaxies(ETGs) ([Taehvun et al.l [2012 ) 
suggest that structural properties such as the effective ra¬ 
dius (radius at which half of the total light of the system 
is emitted) and their surface brightness are quite simi¬ 
lar from those of the other ellipticals or the spheroids of 
galaxies with few or no signs of tidal interaction. Other 
results have shown evidence of ki nematically dec o upled 
components in early-type galaxies (|Hau &; Duncanl [2006; 
iKoorolin &; Zeilingen l20Qol b Given the standard paradigm 
of CDM it is expected that even these “passive” galaxies 
also formed by galaxy mergers, then the similarity of the 
bulge regions for the ETGs could be due to the some mi¬ 
nor mergers (galaxy mergers of unequal mass with mass ra¬ 
tio of 1/10 — 1/100) slightly affecting the bulge but mas¬ 
sive enough to still induce tidal debris, or by major merg- 
ers(mergers of comparable mass galaxies) that happened a 
long time ago so that now the bulge has relaxe d. Cosmolog¬ 
ical N-Body simulations in CDM have shown ([Klypin etahl 
l201ll : iBehroozT et al.l [2oTll l that merger events are frequent 
when galaxies are growing, however, it must be taken into 
account that anticipating the final structure around a galaxy 
is rather difficult given the dependence on the environ¬ 
ment and merger history in the formation of each particular 



Observationally, it is not uncommon to run into difficul¬ 
ties identifying and studying tidal structures around galax¬ 
ies, sometimes these structures are too faint compared to 
the central light distribution that a detailed analysis proves 
impossible. Nevertheless, there are observations where faint 
structures such as shells are clearly seen around galaxies 


at some ] 
ated host( 

particular radii from the center of their associ- 

Turnbull,Bridges V Carter 1999 

•iMalin V Carter 

1983; 

But a et al. 201(|^Laurikainen et al. 

201(1: iFort et al. 

1986; 

Schweizer V Seitzer 20081). In some 

situations it has 


been possible to identify younger stellar popula tions in the 
outskirts compared to those in the inner bulge (| k ang et al.l 
l2009l j probably reminiscent of a late gas accretion. It is worth 
noticing that in the CDM model the spatial location of these 
remnants of past interaction is not predicted and has to be 
tracked with numerical simulations as it depends on the his¬ 
tory of each galaxy. 

Metallicity is another powerful tool to constrain galaxy 
evolution across cosmic epochs as it is known to be cor¬ 
related to gas infall, outflows due to supernovae, and 
star formation rate within the galaxy. Recent studies 
of star form ation focused on local group dwarf galaxies 
(IWeisz et al. 12014a 3) r evealing radial gradients in metal- 
licitv ([Kirbv et all 20131 b and most of them seem to have 
old s tellar populations formed ~ 10 Gyrs ago (|Madau et al.l 
2014 b followed by a more quiescent star formation likely due 
to loss of gas mass at the reionization epoch. There are high 
redshift galaxies (z ~ 3) that also show metallicity gradi¬ 
ents, their regions of high star formation rate(SFR) are also 
of low er metallicitv ([Mannucci et aLll2009l . l20ld : ICresci et al.l 
1201(f ) due to a large g as infall, they a r e sorr ounded by a 
more enriched disc, in iMannucci et 3 <l201(tl the authors 
concluded that these observations are consistent with ac¬ 
cretion of metal-poor gas in massive galaxies where a high 
SFR can be sustained without requiring frequent mergers, 


this pict ure was also suggested by studies of gas rich disks at 
z ~ l-2 (|Forster Schreiber et al.1 (20091 : pTacconi et~ai1l2010l b 
In order to describe the debris around galaxies and 
stellar population it is necessary to include the gas 
dynamics in simulations, Even though there is an uncer¬ 
tainty due to our lack of understanding of the baryonic 
processes, their addition is also required to partially 
reduce some discrepancies of the CDM m odel, for in¬ 
stance, the cusp-core problem(see Ide Blok! (2010) for a 
review), the overabundance of sa.tel l ites(lK KjiLn_eL-aLl 1999; 
Moore et a j] Il999l: ISimon &; Gehal [20071: iBelokurov et all 


20ld : Mack) et al] 120121 : ICarrison-Kimmel et al.l 
or the Too-Big-to - Failfeovlan-Kolchin et al.l 
ICarrison-Kimmel et al.l l2014b f issue. However, 


2014 ) 

20 111 ; 


of these discrepancies might also be solved assuming 


instance, scalar field dark matter (SFDM) (Sin 19941; 

Lee V Kohl 19961: |Guzman & Matos 200(|; 

Matos V Urena 

200lib strongly self-interacting DMfSperg 

el & Steinhardt 

2000; | Vo gelsberger. Zavala V Loebl 

20l3; 

Colin et all 20021: iRocha et al. 2013h. 

warm dark 


It is of particular interest for us the SFDM alternative, 
here the mass of the field is assumed to be very small (~ 
10 -22 eV/c 2 ) such that its de Broglie wavelength is of order 
~ kpc, relevant for galactic scales. The quantum behavior of 
the field has created much interest in the model due to its 
success to account for some discrepancies mentioned above 
with dark matter properties only, for example, the small 
mass keeps the central density from increasing indefinitely 
due to the uncertainty principle in contrast to CD M simulnr 
tions w here super n ova feedba ck is requir ed(Governato et al 


201(1120121: IPontzen fe Govern ato 2012; IScannapieco et al 


2013l b note that the amount of feedback required to 


produce a constant central d ensity may be in co nflict 
with other observations (| Kuzio de Narav V Soekkensl 1 20 111 : 
IPenarrubia et all 1 20121 ; iBovlan-Kolchin et al.l 2014 b 

An interesting property that has been noted for scalar 
field haloes in excited states or in a superposition of 
them is the appearan c e of “ ri pples” in their mass density 
p rofi l esdSeidel &; Suen Il990i: iBalakrishna. Seidel. &; Suen 


1998; Guzman V Urena-Lopez 2004; Urena-Lopez V Bernal 

2Qlol; Bernal et al.l I 2 OI 0 I ; iMatos &; Urena-Lopez 120071 : 

Rob les &; Matosll2013l ) . This property motivated us to inves¬ 
tigate whether features like rings or shells in non-interacting 
galaxies hosted by multistate scalar field haloes could arise 
independently of the galaxy merger mechanism. 

In the next section we give a brief overview of previous 
works in the scalar field dark matter model and describe 
a model for galaxy formation in SFDM haloes. We modify 
the hydrodynamics code ZEUS to implement an evolving 
SFDM halo in which a gaseous component is present, we 
provide the details of our simulations in Section 3. In Section 
4 we discuss our results and section 5 is devoted to our 
conclusions. 


2 SFDM 
2.1 Overview 

The main idea in the scalar field dark matter mode l(Sin 
Il994l : I Ji &; Sinl Il994l : I Lee &; Kohl Il996l : iGuzman &; Matosl 
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12000 : iMatos &; Urenall2001l : iHu. Barkana.fe Gruzinov.11200(1 ) 
considers a self-interacting scalar field with a very small 
mass, typically of ~ 10 _ 22 eV/c 2 , such that the quantum 
mechanical uncertainty principle and the interactions pre¬ 
vent gravitational collapse in self-gravitating structures, 
thus the haloes are characterized with homogeneous den¬ 
sities (usually refered as a cores) in their centers, in gen¬ 
eral the core sizes depend on the values of the mass a nd 
the s elf-interacting param eters (I Colp i. Shapir o. &; Wassermanl 
198fihffor _a review see ISuarez. Robles fc Matosl 11201 .^ : 


Rindler-Daller fc Shapirol (l2014h 'l. 

There has been several studies of the cosmo¬ 
logical evolution that results for a scalar field mass 


evolution of CDM can be reproduced( 

Matos & Ureha 

2001 ; 

Chavanisl 20111: Suarez & Matos 2011 

: iMasana et al. 

2012 ; 

Schive. Chiueh. & Broadhurstl 2014h. 

there is consistency 


with the acous tic peaks of the cosmic microwave background 
radiation (|Rodrfguez-Montova et all l2010h and it implies a 
sharp cut-off in the mass power spectrum for halo masses 
below 1 O 8 M 0 suppr essing structure formation of low mass 
dark matter haloes (iMarsh fe Silkl 120141: Bozek et al. I20l3 : 


IHu. Barkana.fe Gruzinovl l2000l : IMatos fe Urenal ~ 200lh . 


Moreover, there is particular interest in finding equilibrium 
configurations of the system of equations that describe 
the field (Einstein-Klein-Gordon system) and of its weak 
field approximation (Schrodinger-Poisson(SP) system), 
different authors have obtained solutions interpretated 
as boson stars or later as dark matter haloes showing 
agreement with rotation curves in galaxies a nd velocity dis- 


persion promes m a 
199l]: |Lee & Koh 

warr spneroiaai gaiaxiesoeiaei 
19961: Bohmer & Harko 

ouen 

2007; 

Robles & Matosl 

2012, 120131; Lora et al.1 

2012 ; 

Lora & Magana 201^ Martinez-Medina. Robles 1 & Matos 


far the large and small scales observations are well de¬ 
scribed with the small mass and thus has been taken as 
a prefered value but the precise values of the mass and 
self-interaction parameters are still uncertain, we believe 
tighter constraints can come from numerical simulations 
(|Schive. Chiueh. fe Broadhurstl[2014l i and comparisons with 
large galaxy samples which we plan to do in the future. 

Recently the idea of the scalar field has gained 
interest, given the uncertainty in the parameters the 
model has adopted different names in the literature de¬ 
pending on the regime that is under discussion, for 
instance, if the interactions are not present and the 
mass i s ^ 10 ~ 22 eV/c 2 this lim i t was called fuzzy dark 
matter (jhlu. Barkana fc Gruzinov. 2000 ) or mor e rec ently 


wave dark matter ([Schive, Chiueh. &; Broadhurstl i2014b an¬ 
other limit is when the SF self-interactions are described 
with a quartic term in the scalar field potential and dom¬ 
inate over the mass (q u adratic) term, this was studied in 
(|Goodmanl l2000l : ISlepian &; Gooodmanl 12012 ) and called r e- 
pulsive dark matter or fluid dark matter bv( Peebleal2000 i. 


Notice that for a scalar field mass of ~ 10 _ 22 eV/c 2 
the critical temperature of condensation for the field is 
Tcrit ~ m ~ 5 / 3 ~TeV, which is very high, if the tem¬ 
perature of the field is below its critical temperature 
it can form a cosmological Bose Einstein condensate, if 
it condenses it is called Bose-E i nstein condensed (BEG) 
dark matter (|Matos &; Urehall200ll : I Guzman &; Matosl feoOO; 


Berna]_eta 1 2 0inl: iRobles fc Matosl 2013: Harkol 1201 ll : 


Chav anis 201 ll: lLi~ et al.lEoiH l. ISikivie fc Yand (feOOflFl men- 
tioned that axions could also form Bose-Einstein conden¬ 
sates even though their mass is larger than the previous 
preferred value, n o tice t hat the result was contested in 
I Davidson &; Elmier! (2013 ). this suggest that the condensa¬ 
tion process should be study in mor e detail to co nfirm it 
can remain as BEC dark matter. In lUreha- Lonezl ll200flh . 
it was found that complex scalar field with ra< 10 - 14 eV/c 2 
that decoupled being still relativistic will always form a cos¬ 
mological Bose-Einstein condensate described by the ground 
state wavefunction, this does not preclude the existence of 
bosons with higher energy, particulary in dark matter halos. 

We see that the smallness of the boson mass is its char¬ 
acteristic property and cosmological condensation is a likely 
consequence. The preferred mass of the scalar field dark 
matter lies close to ~ 10 _ 22 eV/c 2 satisfying the above con¬ 
straint, although there are still uncertainties on the mass 
parameter, in order to avoid confusion with the known ax- 
ion and help with the identification in future works, we find 
it useful and appropiate to name the scalar field dark mat¬ 
ter candidate, given the above characteristics we can define 
it as a particle with mass m< 10 _ 14 eV/c 2 , being commonly 
described by its wavefunction we choose to name this DM 
candidate psyon. 

It is worth emphasizing that despite the variety of 
names given to the model the main idea mentioned above 
remains the same, it is the quantum properties that arise 
due to the small mass of the boson that characterize and 
distinguishes this paradigm, analoguous to the standard cos¬ 
mological model represented by the CDM paradigm whose 
preferred dark matter candidates are the WIMPs (weakly in¬ 
teracting massive particles), one being the neutralino, we 
see that all the above regimes SFDM, Repulsive DM, Ax- 
ion DM, or any other model assuming an ultra light bosonic 
particle comprise a single class of paradigm, which we call 
Quantum Dark Matter^ QDM) paradigm. As pointed before, 
in the QDM paradigm the small mass of the dark matter 
boson leads to the posibility of forming cosmological con¬ 
densates, even for axions which are non-t hermally pro duce d 
and have masses in 10 -3 — 10 _ 6 eV/c 2 (|Sikivie &; Yanell2009h . 
from here we can obtain a characteristic property that dis¬ 
tinguishes these dark matter candidates from WIMPs or 
neutrinos, namely, the existence of bosons in the condensed 
state , or simply BIOS , from our above discussion the axion 
and psyon are included in the BICS. 

2.2 Dark matter haloes of scalar field 


There has been considerable work to find nu¬ 
merical solutions to the non-interacting SFDM in 
the non-relativistic regime to m odel sphe rically 
symmetric haloes (IGuzman & Urena^Loj^e^ 


Urena-Lopez fc Bernall |201(J: Bernal et al 


Bravl [20121: Ruffini & Bonazzola 196§; iKaupl 


2004; 


2010 ; 


19681 : 


Seidel fe Suenl Il99ll : iLee &; Liml l 2010 |b and a ls o for 
the self-interacting SFDM (iBohm er fe Ha rkol 2007; 


1986; 


Robles.J^lMQj l20l5: Col^kSha^iro^^Wsserman _ 

Rindler-Daller &; Shapiroll2012l : lBalakrishna. Seidel, fe Suenl 

19981 : iGoodmanl l 2000 h. it is w orth notin g that as men¬ 
tioned in IGuzman fe Urena-Loped (|2004h for the weak 
field limit of the system that determines the evolution of 
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a spherically symmetric scalar field, that is, the Einstein 
and Klein-Gordon equations, for a complex and a real 
scalar field the syste m reduc es to the Schrodinger-Poisson 
eauations dArbev et al.1 2 003h . From current bounds re¬ 
ported in iLi et al.l (|2014h obtained by imposing that the 
SF behaves cosmologically as pressureless matter(dust) we 
obtained that the interacting parameter would be extremely 
small for the typical mass of ~10 - 22 eV/c 2 , therefore we 
expect that solutions to the SP system with no interac¬ 
tions would behave qualitatively similar to those when 
self-interactions are included, this assumption is supported 
by the similarity in the sol uti ons for a small self-co upling 
found in other work s (Balakrishna, Seidel, & Suen 1998; 
IColpi.Shapiro.&; Wassermanl 19861 : iBriscesa 201 lh . in this 
work we will study the non-interacting solutions. 

One characteristic feature of stationary solutions of the 
form ^(x, t) — e~ lEnt 0(r) for the SP system is the appear¬ 
ance of nodes in the spatial function <f>(r), these nodes are 
associated to different energy states of the SF, the zero node 
solution corresponds to the ground state, one node to the 
first excited state, and so on. These excited states solutions 
fit rotation curves(RCs) of large galaxies up to the outer¬ 
most measured data and can even reproduce t he wiggles 
seen at large radii in high-resolution observations (Sin 1994; 
IColpi. Shapiro.&; Wassermanl 1 19861 : iRobles &; Matosl l2013h . 
However, haloes that are purely in a single excited state 
seem to be unstable when the number of particles is not con¬ 
served (finite perturbations) an d decay to t he ground state 
with di fferent decay rates (Guzman fe Urena-Lopezl 1 20041 : 


iBalakrishna, Seidel. &; Suen Il998h . though they are stable 
when the number of particles is conserved (infinitesimal per¬ 
turbations). The ground state solution is stable und er finite 
perturbations and infinite simal perturbations ([Bernal et ahl 
l 2010 l : ISeideT &; Suenlll990h . but has difficulties to correctly 
fit the rotation curves in large galaxies because its associ¬ 
ated RC has a fast keplerian behavior shortly after reaching 
its maximum value unable to remain flat enough at large 
radii. 

One way to overcome this problem was to consider 
that bosons are not in one state but instead coexist in dif¬ 
ferent states within the halo, given our intention to de¬ 
scribe dark matter haloes we will refer to such config u¬ 
rations as multistate halos(MSHsl(IUren a-L 6 pe z &; Bernal 
120101 : iMatos &T U rena-L6pezll2007l : IRobles &; Matosll2013l lbh. 
The size of the MSH is determined by the most excited state 
that accurately fits the RC for large radii, excited states are 
distributed to larger radii than the ground state, and in 
contrast to the halo with single state there are MSHs that 
are stable under finite perturbations provided the ground 
state in the final halo confi guration has enoug h m ass to 
stabilize the coex isti ng statedU rena-Lo pez Sz Bernall 120101 : 
iBernal et al.ll2010h . In lBernal et al.l (|2010h it was shown that 
MSHs can be studied in the classical approach as a collection 
of classical scalar fields coupled through gravity, one field ifii 
for each state, this would modify the SP system such that 
its source of energy density woul d be the sum of densities 
in each state (|Urena-Lopez fc BernafeOlO b where each state 
satisfies its respective Schrodinger equation and the states 
are coupled through the Newtonian gravitational potential 
U. 

In spherical symmetry (in units tt= 1, c=l) the SF solu¬ 
tions, <I> n , satisfying the Einstein-Klein-Gordon system are 


related to the SP wavefunctions through 

V87r£M> n (x, t) — e _ * mt '0n(x, t) (1) 

with m the mass of the SF. For stationary solutions we can 
assume wavefunctions of the form ^ n (x, t) = e~ lEnt (f>{r ), 
then the SP system reads 

V 2 0n = 2 (U — E n )4>n 

V 2 t/ = y>n | 2 ( 2 ) 

n 

with V 2 the Laplacian operator in spherical coordinates, and 
E n the energy eigenstates. 

We consider the case of a MSH with only the 
ground and fir s t excited sta tes given that for this MSH 
lUrena-Lopez fc Bernall (12010 ) reported the critical value for 
its stability under small perturbations. If the ground state 
has N ^ particles and there are N ^ in the excited state, 
the authors found that MSH would be stable under small 
perturbations provided the fraction 

jV (2) 

V ]y(iy ^ 1*2 — Vrnax- (3) 

In all the configurations studied in which 77 is larger than the 
maximum for stability rimax , the induced instability causes 
the configuration to lose a few particles and a remarkable 
effect takes place, the excited state does a fast transition 
to the ground state and viceversa, that is, the populations 
of the states are inverted such that the final 77 complies 
with the stability condition ©, after the inversion the MSH 
approaches a stable configuration where the particle num¬ 
ber remains constant for each state, after the transition the 
particles that were initially in the excited state are now in 
the ground state where they have redistributed and become 
more compact than in their original distribution. 

This population inversion(PI) is a characteristic fea¬ 
ture of MSHs and as we will show it can have conse¬ 
quences in the gas distribution and on the accretion rate 
during the period of galaxy formation. In fact, given the 
cut-off in the mass power spectrum due to the small mass 
of the psyon we expect a delay in the first s tructures to col ¬ 
lapse compared to the thos e found in GDM (IMatos &; Urenal 
l200ll : iMarsh &; Silk! 120141 : iBozek et al.l 120141 ). this was 
confirmed with cosm ol ogical simulations for a mass of 
8 x 10~ 23 ey (|Schive, Chiueh. fe Broadhurstll2014h . therefore, 
this intrinsic change of states within MSHs at high redshift 
might have observable consequences in the gas and stellar 
properties of the galaxies that will form in such potential 
wells. 


2.3 Galaxy formation scenario in SFDM haloes 

Initial fluctuations that grow due to the cosmological expan¬ 
sion of the universe eventually separate from it and start 
collapsing due to its own gravity, at this time (known as 
turnaround) the halo has a number of psyons that can be in 
different states, their values would be determined before the 
collapse of the configuration, and have some dependence on 
the its local environment. Depending on the number density 
of bosons populating the excited states we can have different 
fates for the haloes as mentioned in subsection ( 2 . 2 ). 

From the mentioned results of the literature and based 
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on the successful fits to different galaxies we propose the 
following model of galaxy formation in SFDM haloes. 
The smallest and less dense systems, such as dwarf galax¬ 
ies, would reside in SFDM haloes with most bosons in 
the groun d state except possib ly for just a few excited 
particles (|Martinez-Medina. Robles. &; Mat osl 1201 5h . this is 
because the potential wells of these haloes, being the lowest 
dense systems, would be just massive enough to collapse 
and allow the existence of the state of minimum energy that 
would form a bound configuration. Larger configurations 
that had initially a larger number of bosons in excited states 
than in the ground state can undergo a population inversion 
and reach a stable state, collapse to a dense ground state 
or become a black hole depending on how large the fraction 
of particles in excited states is after turnaround, given the 
different possible outcomes for this case we expect that 
most galaxies are formed this way, for instance, ellipticals 
can form when haloes transition to a stable configuration to 
form a ground state and quickly deepen the gravitational 
potential resulting in a dense and compact structure where 
the bulge can form, for high surface brightness galaxies 
which usually have RCs that fall slowly after its maximum, 
these galaxies would likely correspond to SFDM haloes 
whose final stable configuration has a comparable number 
of bosons in the excited and ground states as both states 
would contribute evenly in the central region increasing 
the potential well where more star formation can take 
place, and for distances larger than a given radius (likely 
the first maximum peak in the RC) the ground state 
would be a small contribution to the density, thus in the 
outer region the baryonic RC would remain below its 
maximum. Another possible outcome for extended haloes is 
when the fraction of excited particles is below its stability 
threshold, here gravity will slowly redistribute DM bosons 
in the stable MSH but mostly keeping the particle number 
constant, considering that psyons in the excited states are 
subdominant but more widely distributed than the ones in 
the ground state of a MSH, the DM gravitational potential 
in the center would be less dense than in the high surface 
brightness case so that gas accretion would also proceed 
more slowly resulting in rotation curves for the baryonic 
matter that slowly increase to large radii or with almost 
flat profiles, this case is quite simila r to what is fo und 
low surface bri ghness galaxies feuzio de Narav et al.l 


20ld: 


Kuzio de Narav fc Spekkend 


Kuzio de Narav fc Kaufma.nril 201 il l . 


HH; 


This galaxy formation scenario broadly describes the 
different galaxy types showing that only due to the dark 
matter properties in the context of SFDM it is possible to 
agree with the general features of several galaxies. However, 
we do not expect this scenario to represent all galaxies, in 
fact, to get a better description of individual galaxies we 
require taking into account other fundamental parameters 
that affect their evolution, such as the environment, angular 
momentum, galaxy mergers etc., we will leave a more in dept 
exploration of the scenario including the full astrophysical 
processes for a future work. 

For now, we take the given description of elliptical 
galaxies to explore whether some substructure found in some 
of these galaxies can be a consequence of the properties of 
the background DM halo where such ellipticals are formed. 
We do not discard the possibility that shells or rings result 


- transition of DM states occurs at lGyr 


Time evolution of p gas -1 


of DM states occurs at lGyr 



Figure 1 . Density profiles for dark matter and gas for simu¬ 
lations A \ (upper panels) and Af (lower panels). Shown are the 
initial(cyan lines) and final(yellow lines) DM profiles of the MSH 
composed of the ground and first states of the scalar field. Also 
shown are the initial condition for the disc of gas (gray lines), 
its density and rotation curve just before (blue-dashed line) and 
shortly after (red-dashed line) the population inversion occurs in 
the SFDM halo to reach its stable configuration, the black dashed 
line shows the profiles at the end of the simulations when the mul¬ 
tistate halo is now stable. 


from mergers, instead, we explore a new possible mechanism 
in which they can arise offering an alternative explanation 
to the tidal debris in certain non-interacting ellipticals. We 
study this issue in a multistate halo. 


3 SIMULATIONS 

We consider a multistate halo with only the first and ground 
state coexisting, this is because we know the stab ility thresh¬ 
old for this multistate configuration. In lBernal et al.l (2010) 
it was seen that the larger the r) above the critical value, 
the faster it settles to its final stable configuration, to be 
sure that the MSH is initially in the unstable regime we 
pick an intermediate value rj— 1.6. We explore the evolu¬ 
tion of the gas distribut ion that follows a disc of g as with a 
Miyamoto-Nagai profile ([Miyamoto &; Nagailll975l i with gas 
mass of M g m 3.9 X 10 9 M©, and for the horizontal and 
vertical scale-lengths we use a — 6.5 kpc and b — 0.5 kpc 
respectively (we also simulated some cases with the same 
gas mass but using an initial spherical Plummer profile but 
found no change in the main results, there was only a slight 
increase in the computation time, we then report only the 
cases for an initial disc of gas). 

The code we use to evolve the gaseous component 
emb edded in the MSH is the hydrodynamics code ZEUS- 
MP (|Haves et al.l[2006 ). It is a fixed-grid time-explicit Eule- 
rian code, we solve for the standard hydrodyn amics e qu ation 
for the baryonic component as in lMartinez-Medina Matosl 
ll20l4l . but here, we modify the code to allow for the evo¬ 
lution of a multistate dark halo that undergoes a transi¬ 
tion as described below. In order to implent the popula¬ 
tion inversion in ZEUS we use a semi-analytical model that 
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within Rj is 

M A R A=[l) P -^- (7) 

During the population inversion in the numerical so¬ 
lution both st ates lose a few percen t age o f their initial 
masses(<10% ([Urena-Lopez &; Bernall l20ld : iBernal et al.l 
2010)), but this loss would not substantially affect the final 
halo profilea, for this reason we will neglect the mass loss 
in our analytical study of the MSH. In order to construct 
the initial configuration in an unstable state we work with 
a mass ratio r\ — 1.6, when the population inversion occurs 
the initial mass ratio M^R^)/M{(R\) will also be inverted 
and the final MSH will be stable. After the inversion the 
total initial mass of the ground state M{(R\ ) becomes the 
final mass of the excited state M^R^) and viceversa, and 
the nodes in the wavefunctions will also have changed, that 
is 


Figure 2 . Density profiles for dark matter and gas for simu¬ 
lations A\ (upper panels) and A\ (lower panels). Shown are the 
initial(cyan lines) and final(yellow lines) DM profiles of the MSH 
composed of the ground and first states of the scalar field. Also 
shown are the initial condition for the disc of gas (gray lines), 
its density and rotation curve just before (blue-dashed line) and 
shortly after (red-dashed line) the population inversion occurs in 
the SFDM halo to reach its stable configuration, the black dashed 
line shows the profiles at the end of the simulations when the mul¬ 
tistate halo is now stable. 


let us control the rate of the transition to a stable state 
helping with the interpretation of the results. To achieve 
this we first use an analytical approximation to the nu¬ 
merical solutions obtained in the Newtonian limit for a 
confined distri bution where the gravitatio n al potential is 

weak (I U rena-Lopez fc B ernafeOldnGuzman fe Urena-Lope3 

12004 iRobles &; MatosH2013h . the scalar field density profile 
due to the wavefuntion in the state j is 


PjiT) = p c ,j 


sin 2 (inr / RA 
(jirr/Rj) 2 


(4) 


where j—1+n labels the state of the scalar field and 
n=0,1, 2,... is the number of nodes in the profile, thus j — 1 
denotes the ground state with zero nodes and j — 2 cor¬ 
responds to the first excited state, and the total density of 
the MHS would be described by p(r)=p±(r)-\-p2(r) . In equa¬ 
tion © Pc, j represents the central density of state j and Rj 
is taken as a truncation radius of the corresponding state, 
that is, for radii grater than Rj the number of particles in 
the state j is neglected, thus we take pj(r)= 0 for all r > Rj 
whereas for r ^ Rj the density profile of the state j is given 
by (gj). 

The mass profile and rotation curve for state ,7 are 
([Robles &; Matosll2013l : iMartinez-Medina fe Matosll2014| j 


M,(r) 


v f M 


4tt pcj r ( _ sin(2fcj-r) \ 

k) 2 y 2 kjr J ’ 

47 rGpcjf^ sin(2 kjr)\ 

2 k 2 \ 2kjr J 7 


(5) 

(6) 


where we defined kj := jir/Rj. The total mass enclosed 


Mi(Ri) —► Mf(Ri), Mi(Ri) —»■ M[(R{) (8) 


j: 1,2 —>2,1=/. (9) 

We still need to specify a function to model the time evo¬ 
lution of the transition to the stable configuration. As we 
mentioned, for unstable haloes the transition happens re¬ 
ally fast, during this period the central density varies until 
it reaches the new stable value, the inversion also modifies 
the radius of each state. In our semi-analytical approach the 
radius and density are related to the mass through eq. ©, 
thus we can choose to vary the radius to induce the transi¬ 
tion, remembering that the transition is relatively fast, we 
found that a smooth function that captures the PI of the 
numerical evolution and allow an interpolation between an 
initial and final Rj is 


Rj, (t) = R) + ( 1 + tanh[a(t - t inv )\ 


( 10 ) 


where a determines the transfer rate from the initially unsta¬ 
ble to the final stable MSH, €j,j = R*., — Rj is a parameter 
that relates the initial radius to the one after the inversion, 
and ti nv is the time where the PI is halfway to reach its fi¬ 
nal stage with j' confined in R(, , we applied a corresponding 
function to invert j to j'. 

Using this semi-analytical inversion model we can iden¬ 
tify and study the dependence of the resulting properties 
of the gaseous component on the transition time, in par¬ 
ticular the effect on the production of tidal features. We 
modify the code ZEUS to implement the model and run 
simulations for a Miyamoto-Nagai disc of gas, the gas has 
no self-gravity and is used as a tracer of the gravitational 
potential, the initial condition for the gas velocities are set 
directly from the background dark matter halo and there¬ 
fore match the i nitial halo velocities(see Fig. 1) . As noted by 
other au thors(| Guzman fe Urena-Lopez 120041; lUrena-Looea 


120021 : lUrena-Lopez.Valdez-Alvarado. &; Becerrill 2012h sta- 
tionary solutions to the SP system (for a specified number 
of nodes) are related to each other by a scaling transforma¬ 
tion, so the properties found here for our simulated parame¬ 
ters are expected to appear in similar configurations under a 
suitable scaling parameter. Table [I] shows the parameters of 
mass and radii for our MSH before and after the transition, 
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Table 1. Initial and final parameters defining the MSH. 



Initial 


Final 



3 


r 



1 

2 

1 

2 


1.844 

2.953 

2.953 

1.844 

R * 

10 

20 

10 

20 


a Mass in units 1 X 1O 1O M0 
^ Units in kpc 

these are the values we will use that define the background 
MSH in all our runs, we explore an early and late transition 
specified by ti nv =l and 3 Gyrs respectively, and we also 
evolve these cases when a = 1,1/4 corresponding to an fast 
and slow inversion, we denote by A a realization that 
had a PI with parameters ti nv and a. 


4 DISCUSSION 

In Figure 1 we show the density profiles and the rotation 
curves for the MSH and the disc of gas comparing the cases 
of a rapid early and late transitions (A\ and A\). We notice 
the decline in the gas RC in both cases, in fact, the behavior 
is quite similar in both simulations, the gas moves according 
to the potential of the MSH and when the PI takes place it 
induces a rapid gas infall towards the center (see red dashed 
line in Fig. 1) seen as an bulge like region of high gas density. 
In Figure 2 we plot a zoom to the gas densities of Figure 1 
where we observe a conspicuous peak that appears as a result 
of the halo transition, soon after the transition a portion of 
the outer gas is pulled to the center and gets redistributed 
again in the now stable MSH configuration but now some of 
the gas remains trapped within the region of the first node of 
the MSH density profile, this results in the notable spike in 
the gas profiles. This density increase appears in both early 
and late fast transitions showing that this feature is expected 
independently of when the transition to the stable multistate 
halo happened. We verified that this is not a resolution effect 
by simulating A\ for three different grid resolutions giving 
the density profiles shown in the bottom panel of Fig. 2, the 
inner profile does change due to the higher resolution that 
captures more details of the gas while the spike at ~ 9 kpc is 
only slightly modified and at the location of the node in the 
MSH, suggesting that its origin lies in the ripple-like feature 
of the background MSH. 

This surprising agglomeration of gas at 9 kpc is nev¬ 
ertheless small compared to the mean central density(~ 
0.45M©pc -3 ) being ~ 20% at the peak, the fact that it re¬ 
mains well after the transition is important for observations 
as this can be the origin of some structures around galax¬ 
ies. In fact, in Fig. 3 we plot the gas spatial distribution at 
the beginning and at the end of the simulation with early 
and late halo transition, that is A\ and Af, and we read¬ 
ily observe the presence of a sorrouding ring of gas orbiting 
around the denser central concentration reminiscent of those 
found in certain elliptical galaxies. The ring keeps its loca¬ 
tion until the end of the simulations in our isolated galaxy 


because its origin lies in the shape of the background SFDM 
halo, moreover, the stability of the halo under small pertur¬ 
bations, such as accretion of small galaxies, mildy change 
the inner DM potential but may alter the ring structure 
in a similar way that the population inversion modified the 
initial gas profile, however, once the merging process has 
concluded the gas will start redistributing and according to 
our results had enough gas remain on the outskirts a ring¬ 
like structure would appear again, the final amplitude of 
the ring’s density will depend more strongly on the details 
of the merging process but its location should be similar to 
the non perturbed halo provided the DM halo potential was 
not severely modified, as expected for minor mergers in the 
past or those relatively more recent, thus, tidal features are 
not exclusive of galaxy mergers when the dark matter is an 
ultra light scalar field. 

Comparing the early and late transitions in the MSH we 
obtain similar structures, thus we may focus on one of them, 
we choose the case of a late population inversion(U n ^=3 
Gyrs) to explore the effect on the ring’s shape under a 
fast(<a=l) and slow (a=l/4) transition. Fig. 4 depicts the 
comparison at the end of the simulation, we observe that 
both Af and A\/ a form the ring as expected but the slow 
transition spreads the gas around the ring creating a broader 
substructure, this is consistent with our galaxy formation 
scenario, as the gas accretion proceeds slower for A\^ the 
ring will slowly form from the concentration of consecutively 
infalling gas, on the other hand in A\ the rapid transition 
enhances the gas infall and accelerates the gathering of gas 
at the ring’s location quickly depleting the rest of the gas in 
low density regions creating a thin and more confined ring. 

Moreover, a fast gas concentration is likely to induce a 
period of star formation, a massive accretion in the early uni¬ 
verse would generate a starburst producing several stars and 
winds that could drive gas away from the center, part of the 
latter will be recycled and some can be trapped in the loca¬ 
tion where the outer ring forms, there the gas can trigger star 
formation too or just accumulate as dust, the central star 
formation affects the metallicity in that region and if the gas 
infall stops early due to effective gas blowouts consequence 
of a high star formation rate the galaxy would be left with a 
dominant fraction of old stars. For a massive galaxy, part of 
the ejected gas cannot escape and will remain bounded scat¬ 
tering the inner photons and producing a luminous bulge, 
this is consistent with the formation of some ellipticals, in 
particular this picture could describe the puzzling origin of 
dust rings in the Sombrero galaxvfM104j( Bai aia et al.ll l984; 
lEmsellemlll995h . although a more careful study of this issue 
has to be done. 

The gas fraction and nature of the dark matter halo 
are fundamental to determine the fate of a galaxy, for in¬ 
stance, dwarf galaxies reside in SFDM haloes in the ground 
state and have shallow potentials, then if a starburst occurs 
early it can blow out a large portion of gas that if it is ac¬ 
creted much later the next generations of stars will differ 
from the old population observed as two different fractions 
of stellar population where the younger one is located at 
larger radii due to its late accretion, analyses of stellar age 
and metallicity at different radii can be a way to distinguish 
between CDM predictions and the SFDM, future hydrody¬ 
namics simulations of SFDM will be able to confirm these 
predictions. 










8 Victor H. Robles, L. A. Martinez-Medina and T. Matos 



Figure 3. Comparison of initial gas distribution(left) and its final stage. The central panel corresponds to a halo that had an early 
transition to its stable configuration at ti nv =1 Gyr(A^) and on the last panel(right) the transition was at a later epoch ti nv = 3 Gyrs(A^). 
A ring appears in both simulations at the corresponding density peak of Fig. 2. 


The existence of the ring is a generic feature in our 
simulations of isolated galaxies, notice that for an initially 
stable halo (77 < rjmax) the ring will also appear, it suffices 
to observe that our final MSHs fall precisely in this regime. 
The SFDM model provides a new mechanism to form rings 
without major changes to the galaxy inner region mainly 
because it is where the ground state is mostly confined, in 
fact, the latter can exhibit a solitonic behaviour keeping its 
profile even after several mergers as sh own in SFDM cos¬ 
mological simulations (|S chive et al. 2014b). Although perfect 
rings are created in our simulations, we note that this struc¬ 
ture can be deformed by the passage of an accreted satellite 
or due to tidal forces from nearby galaxies culminating in 
incomplete or displaced rings, t hese tidal debris would man¬ 
ifest in the more common shells(Malin &; Cart er 1980, 1983; 


iTurnbull.Bridges fe Carterlll99Sf 


Galaxies that undergo major mergers form under vio¬ 
lent conditions, thus the shape of their tidal structures can 
vary widely especially when the interaction is still ongoing, 
for these latter galaxies far from equilibrium the emergence 
of an exotic tidal feature might be more accurately explained 
by modeling the merging galaxy and evolving it in a simu¬ 
lation that includes baryonic processes and a live halo. For 
the more isolated systems like the considered here, we have 
found that there is another mechanism that leads to outer 
structures but with the difference that they posses a certain 
degree of symmetry associated to the background DM halo 
and not to a particular merger event. 


It is interesting that the quantum properties of the 
psyons manifest on a macroscopic scale and have implica¬ 
tions in observable quantities like gas and stars, remark¬ 
ably a similar situation where quantum effects manifest 
on larger scales has already been observed in local exper¬ 
iments, in particular it was seen that water droplets can 
bounce indefinitely on the surface of a vibrating fluid bath 
and propel themselves across the surface of the fluid by 
virtue of their pilot-wave dynamics exhibiting quantum fea¬ 
tures like single-particle diffracti on, tunneling , quantized or¬ 
bits, and orbital level split t ing (I Walker 19781: Couder et al. 
2005; I Couder &; Fortl 120061 : IProtiere. Boudaoud. &; Coudei 



Figure 4. Gas distribution after 5 Gyrs for simulations 
A\ (upper-left and its colored scale in the bottom-left panel) 
and A^j 4 (upper-right and its colored version in the bottom-right 
panel). The ring on the left panels is more defined as the MSH un¬ 
dergoes a rapid transition to stability while for a slow transition 
(right panels) more particles are spread around the ring. 


2006 : lEddi et al.l 120091 120121: iFort et alJ l2010l: iBushl l20ld : 
Ozaa, Rosales, & Bush 201-i)) . In the SFDM context the ana¬ 


logue fluid would be the dark matter and the baryons would 
be the droplets guided by the DM, whether this analogy is 
valid deserves further consideration and will be treated in 
other work. 





































5 CONCLUSIONS 

In this work we have seen that if the dark matter is a ultra 
light scalar field the internal structure of galaxies can reflect 
fossil records of their formation, in our case we study the 
formation of tidal structures similar to rings and shells in 
non-interacting galaxies. We discussed the nature of these 
structures in the context of SFDM and conclude that inde¬ 
pendently of galaxy mergers they can arise as a consequence 
of the quantum nature of the scalar field dark matter bo¬ 
son (hereafter psyons) that form the halo where galaxies are 
embedded in, mainly, due to the intrinsic ripples in the dark 
matter halo density profile. 

The possibility of describing tidal substructure in scalar 
field dark matter halos was also suggested in iBravl ( 2012 ) 
but under different halo configurations and without includ¬ 
ing the evolution of a gaseous component, here we make 
use of spherically symmetric haloes in multistates because 
their stab ility has been studied and confirmed in prey i- 
ous works ([ Urena-Lopez fe Bernalll2010l : iBernal et al.ll2010h . 
In addition, they possess desirable properties to describe 
a variety of galaxies. Ins pired by the numerical solutions 
in lUrena-Lopez fe Berna] ( 20 Iff ), where they reported the 
maximum fraction of psyons in excited states (eq.[3]) to form 
a stable multistate halo composed of the ground state and 
first excited state, we propose a semi-analytical model to de¬ 
scribe a similar halo configuration and study the distribution 
of a gaseous disc that is contained within such scalar field 
DM halo. The halo is set to be initially unstable by choos¬ 
ing a fraction of excited boson above the stability limit, this 
induces a transition to reach a final stable configuration af¬ 
fecting the initial gas distribution, after the transition the 
population of the states is inverted and the halo stabilizes 
once the ground state dominates the halo mass that by grav¬ 
ity will be confined later in the center of the halo. At the end 
of the simulations (5 Gyrs) part of the gas keeps orbiting in 
the form of a ring outside the main bulge-like region, the 
structure of the ring is not significantly different when the 
halo transition happens early (ti nv =1 Gyrs) or late(U nv =3 
Gyrs) but the ring broadens if the transition is slow and be¬ 
comes more defined for fast transitions. Due to the lack of 
stability parameter for MSH in higher energy states we ex¬ 
plore the first non trivial superposition of multistates, given 
the positive results it seems that tidal structures would also 
be found at larger radii for SFDM haloes that have boson in 
higher excited states, although they would become fainter 
with distance, this would agree with shells that appear to 
be distribited at larger radii. 

The tidal features in this work have a distinct origin 
from the usual tidal stripping of galaxies, they emerge from 
the quantum properties of the psyons that compose the dark 
halo, the presence of the symmetric ring is due to the spher¬ 
ically symmetry hypothesis, deviations from this symmetry 
could result in incomplete rings and shells although a cer¬ 
tain degree of symmetry should remain which distinguishes 
to remnants that come from recent galaxy mergers that de¬ 
pend on the path of the accreted satellites and does not 
necessarily leave tidal debris of a particular symmetry, our 
results suggest that outer regions of dark haloes may not be 
totally smooth. Exploring in more detail the regions where 
shells and rings form could also help to put constraints on 
their origin, their formation, and serve to test galaxy and 
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halo formation scenarios such as the quantum dark mat- 
ter(QDM) paradigm. 

Additionally, we have provided a classification for the 
dark matter models that are based on the assumption that 
the dark matter is a scalar field of a very small mass that 
proves useful for future references and that helps to es¬ 
tablish the main feature of the particle candidates of this 
paradigm. We formulate a galaxy formation scenario based 
on known results of halo formation in the SFDM context 
that can explain qualitatively the general properties of the 
different galaxy types in the universe and applied to inter¬ 
pret our findings. Given the success of the SFDM model to 
account for the apparent discrepancies of the standard CDM 
model with only dark matter properties, it will be desirable 
to include astrophysical processes in simulations to get a 
fair comparison with state-of-the-art CDM simulations. The 
QDM paradigm is indeed an interesting possibility that de¬ 
serves further study. 
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